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ABSTRACT 

We present near-infrared spectroscopy for a complete sample of 33 ultraluminous 
infrared galaxies at a resolution of R ~ 1000. Most of the wavelength range from 1.80- 
2.20 /um in the rest frame is covered, including the Paa and Br7 hydrogen recombination 
lines, and the molecular hydrogen vibration-rotation 1-0 S(l) and S(3) lines. Other 
species, such as He I, [Fe II], and [Si VI] appear in the spectra as well, in addition 
to a number of weaker molecular hydrogen lines. Nuclear extractions for each of the 
individual galaxies are presented here, along with spectra of secondary nuclei, where 
available. The Paa emission is seen to be highly concentrated on the nuclei, typically 
with very little emision extending beyond a radius of 1 kpc. 

This survey was carried out primarily to search for signatures of active nuclei via 
velocity-broadened hydrogen recombination or the presence of the [Si VI] coronal line. 
These signatures are rare in the present sample, occurring in only two of the 33 galaxies. 
The extinction to the hydrogen recombination lines is investigated via the Paa /Br7 line 
ratio. It is found that visual extinctions to the nuclei in excess of 10 magnitudes are 
relatively common among ULIRGs, and that visual extinctions greater than 25 mag are 
necessary to conceal a QSO emitting half the total bolometric luminosity. The ionized 
hydrogen regions in many ULIRGs are sufficiently obscured that dust-enshrouded active 
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galactic nuclei would remain hidden at 2 at the current level of sensitivity. The 
vibration-rotation lines of molecular hydrogen appear to be predominantly thermal 
in origin, with effective temperatures generally around 2200 K. The relative nuclear 
velocities between double nucleus ULIRGs are investigated, through which it is inferred 
that the maximum deprojected velocity difference is ~ 200 km s _1 . This figure is lower 
than the velocities predicted by physical models of strong interactions/mergers of large, 
gas-rich galaxies. 

Subject headings: galaxies: infrared — galaxies: starburst — galaxies: active 

1. Introduction 

Understanding the nature of the ultimate power source in ultraluminous infrared galaxies 
(ULIRGs) has been a primary driver for numerous scientific investigations ever since the discovery 
in the 1980's that ULIRGs represent the most luminous class of galaxies in the local universe. 
With total power outputs comparable to those of quasars, characterizing the source of power is 
of fundamental importance in understanding the composition and evolution of normal galaxies. 
ULIRGs are found to share the common trait that virtually all have recently or are currently 
experiencing major encounters with other galaxies (Murphy et al. 1996; Clements et al. 1996; 
Sanders et al. 1988; Carico et al. 1990). Encounters like these may play a significant role in 
the process of building normal galaxies. With the discovery of faint sources at sub-millimeter 
wavelengths having characteristic spectral energy distributions like that typically found in ULIRGs 
(e.g., Arp 220) and redshifts z > 2, the importance of understanding these systems and their role 
in the evolution of normal galaxies has been heightened. 

The present survey was designed to look for near-infrared spectral signatures of obscured active 
galactic nuclei (AGN) in ULIRGs, either by the presence of velocity-broadened Paa emission or by 
the appearance of the high excitation [Si VI] coronal line. A summary of the conclusions of this 
survey, along with a median ULIRG near-infrared spectrum, are presented in Murphy et al. (1999). 
Only two of the 33 sample galaxies show clear indications for AGN activity, with the rest appearing 
to be dominated by star formation processes. Either the AGN phenomenon is relatively rare among 
ULIRGs, or the extinction at 2 /jm is still too high to permit a proper assessment. Mid infrared 
studies (e.g., Genzel et al. 1998; Rigopoulou et al. 1999) favor the former conclusion, though it 
is found in this study that the extinction at 2 /im may often be greater than four magnitudes, 
resulting in over one magnitude of extinction at mid-infrared wavelengths. 

This paper presents the individual nuclear spectra for the sample galaxies, for both the pri- 
mary and secondary nuclei where possible. Two-dimensional spectra are also presented to give a 
qualitative view of the line emission distributions typically found in ULIRGs. Some discussion of 
extinction measures, molecular hydrogen excitation mechanisms, and a summary of typical line ra- 
tios found in ULIRGs is included, as is a discussion of the velocity differences found among double 
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nucleus systems. 

2. The Sample 

The sample of 33 ULIRGs in this survey originates from the 2Jy sample of IRAS (Infrared 
Astronomical Satellite) galaxies defined by Strauss et al. (1990, 1992). The galaxies in the present 
sample are chosen to meet the following criteria: 

1. 60 fim flux density, -F„(60 fxm) > 1.94 Jy 

2. Moderate or high quality IRAS flux measurements at 60 and 100 fim (as defined in the IRAS 
PSC, Explanatory Supp. 1988) 

3. F2(60 H> f v( 12 ^m)xi^(25 jum), where F U (X) is the IRAS flux density 

4. Far- infrared luminosity Lfir > 5 x 10 11 L Q 

5. Absolute Galactic latitude, \b\ > 5° 

6. Declination, 5 > —35° 

7. Redshift, 0.055 < z < 0.108 

The first six conditions are those used to define the northern 2 Jy ULIRG sample in Murphy et al. 
(1996), such that the current sample is a subset of this work. The seventh condition places both the 
Paa and Br7 lines in the K band atmospheric window. The far-infrared luminosity is computed 
from the IRAS flux densities by the following conversion: 

Lfir = 3.86 x 10 5 D 2 [2.58-F„(60/im) + F„(100/«n)] L Q , 

(IRAS PSC, Explanatory Supp. 1988, p. X-13) where D represents the luminosity distance to the 
source in Mpc. 

Due to the redshift constraint, the current sample approximately represents a complete volume- 
limited sample, because even the galaxies with the minimum 60 fim flux density are luminous 
enough to meet the sample criteria at all allowed redshifts. At the outer redshift limit of the 
survey, corresponding to a luminosity distance of 455 Mpc using Hq = 75 km s _1 Mpc -1 and 
qo = (as is assumed throughout), the minimum 60 /im flux density of 1.94 Jy by itself accounts 
for 4 x 10 11 L Q of the total far-infrared luminosity. To make up the remaining 10 11 L Q necessary for 
sample inclusion, the 100 /im flux density needs only be greater than 1.25 Jy, implying a maximum 
F u (60 ^m)/F I/ (100 fim) ratio of 1.55. ULIRGs very rarely have 60 /im to 100 ji/m flux density ratios 
this high, such that the volume limit imposed by the upper redshift cutoff is more stringent than 
the flux density cutoff, resulting in a sample that is mostly volume-limited rather than flux-limited. 
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Taking the 64 ULIRGs defined in the northern 2 Jy sample of Murphy et al. (1996), one finds an 
average -Fj,(60 / um)/i ? !/ (100 fim) ratio of 0.90 ±0.22, with a median value of 0.86 and a maximum of 
1.64 (in IRAS 08572+3915), which is the only ratio exceeding the maximum of 1.55 compatible with 
the volume-limited case. The ULIRG from the northern 2 Jy sample with the second highest ratio 
is IRAS 05246+0103, with i^(60 fj,m)/F v (100 fim)= 1.31. Both of these galaxies are contained in 
the present sample. The redshift limit excludes many of the highest luminosity ULIRGs, as these 
are most likely found at larger distances. For instance, Mrk 1014 — a quasar and a ULIRG — would 
be a member of the sample if not for the upper redshift limit. 

Of the 35 ULIRGs in the redshift survey of Strauss et al. (1992) satisfying the above con- 
ditions, two are not included here for the following reasons. The spectra of IRAS 21396+3623 
show this galaxy to be at a redshift of z = 0.1493, rather than the previously reported z = 0.097. 
IRAS 19297—0406 was free of any discernible lines, which may be the result of improper pointing, 
producing the spectrum of a nearby star in the very crowded field. Therefore the entire sample 
contains 33 galaxies, and is very nearly complete in the volume-limited sense. A list of the sample 
galaxies appears in Table 1. 

3. Observations and Data Reduction 

Observations were conducted using the Palomar longslit infrared spectrograph (Larkin et al. 
1996), operating on the 200-inch Hale Telescope and employing a 256x256 HgCdTe (NICMOS 3) 
array. Observations were made over a time period beginning in 1995 July, and ending in 1997 
December. Each object was observed in three grating settings, with each spanning ~ 0.12 fim at a 
scale of ~ 0.0006 fim pixel -1 . Typical slit widths of ~3.75 pixels, or / .'625, resulted in spectra with 
resolutions of R = A/AA « 1000, corresponding to ~300 km s _1 . The grating settings were chosen 
to cover the Paa line at 1.8751 jum, the suite of Br8, H2 1-0 S(3), and [Si VI] lines centered at 1.954 
fim, and the H2 1-0 S(l) and Br7 lines centered at 2.14 fim. The slit was generally rotated such 
that spectra of both primary and secondary nuclei were obtained simultaneously when possible, or 
otherwise to coincide with the major axis of the galaxy, where evident. 

Observational details are provided for each source in Table 2. Individual exposure times 
were 300 s, with total integration times generally around 1800 s. In all cases, the galaxy was 
dithered between two positions on the slit, generally ~ 20" apart, with a smaller scale dither 
pattern employed about these points to eliminate the effects of static bad pixels. In this way, the 
sky integrations were obtained simultaneously, with the effective sky position alternating by 20" to 
either side. Wavelength calibration was provided either by the atmospheric OH airglow spectrum, 
or by arc lamp spectra for the H2 1-0 S(l)+Br7 grating settings, where the OH lines become 
unavailable. Wavelengths reported here are in air units, with the OH wavelength data coming 
from Oliva & Origlia (1992). Atmospheric transmission variations as a function of wavelength 
were compensated via observations of the nearly featureless spectra of G dwarf stars at similar 
airmass, and similar telescope pointing, when possible. The G star observations were performed 
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Table 1. Spectroscopic Survey Sample 



Galaxy Name cz Lib, Separation a 







(km s ) 


(10 L ) 


(arcsec) 


IRAS 


00262+4251 


29205 


12.08 


s 


IRAS 


01521+5224 


23959 


11.95 


5.5 


IRAS 


04232+1436 


23640 


11.99 


s 


IRAS 


05246+01 03 b 


29105 


12.05 


6.0 


IRAS 


08030+5243 


24946 


11.97 


s 


IRAS 


08311-2459 b 


30167 


12.40 


s 


IRAS 


08344+5105 


28999 


11.94 


s 


IRAS 


08572+3915 b ' c 


17491 


12.08 


5.6 


IRAS 


09061—1248 


22014 


11.97 


5.6 


IRAS 


09111—1007 


16247 


11.98 


(37) 


IRAS 


09583+4714 


25748 


11.98 


25.5 


IRAS 


10035+4852 


19449 


11.93 


9.8 


IRAS 


10190+1322 


22898 


11.98 


4.0 


IRAS 


10494+4424 


27671 


12.15 


s 


IRAS 


11095—0238 


31936 


12.20 


s 


IRAS 


1211 2+0305 c 


22009 

Z. _ V / V - / 


12.27 


2.9 


IRAS 


1434g_1447 c 


24798 


12.28 


3.3 


IRAS 


14352—1954 


26938 


11.95 


s 


IRAS 


14394+5332 


31401 


12.04 


(28) 


IRAS 


15245+1019 


22629 


11.96 


2.6 


IRAS 


15250+3609° 


16515 


11.99 


s 


IRAS 


15462-0450 


30060 


12.16 


s 


IRAS 


16487+5447 


31110 


12.12 


3.1 


IRAS 


17028+5817 


31805 


12.11 


13.0 


IRAS 


18470+3233 


23520 


12.02 


(7.0) 


IRAS 


19458+0944 


29983 


12.36 


s 


IRAS 


20046-0623 


25286 


12.02 


s 


IRAS 


20087-0308 


31613 


12.39 


s 


IRAS 


20414-1651 


26061 


12.19 


s 


IRAS 


21504-0628 


23343 


11.94 


s 


IRAS 


22491-1808 c 


23264 


12.12 


1.6 


IRAS 


23327+2913 


32225 


12.03 


12.5 


IRAS 


23365+3604 


19305 


12.13 


s 



a Nuclear separation is given where applicable. S denotes 
a single nucleus. Numbers in parentheses denote a double 
nucleus system for which only the primary nucleus spectrum 
is presented, generally because the separation exceeds the slit 
length. 

b Galaxies with "warm" infrared colors: F^(25 /xm)/i ? „(60 
/im)> 0.2. 

c Contained in the Bright Galaxy Sample (BGS) of Soifer 
et al. (1987). 
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either immediately prior to or following each galaxy observation. Stellar types and airmasses of 
the calibrators are listed for each object in Table 2. The light from the calibrator star was made 
to uniformly fill the slit aperture by chopping the telescope secondary mirror in a triangle-wave 
pattern. As such, the atmospheric calibrators perform the dual function of operating as the flat 
field for the spectral observations. Details about the treatment of the spectral calibration appear 
below. The calibration star observation for IRAS 12112+0305 was compromised beyond repair, so 
the spectra for this galaxy were divided by an atmospheric template. 

Each spectrum is reduced pair-wise in the following manner. Pairs containing spectra at the 
two opposite slit positions are subtracted, with cosmic rays and static bad pixels interpolated. 
The data are divided by the G star spectrum and multiplied by a Planck function, described in 
more detail below. A wavelength reference image is produced using either a median combination 
of the data frames to produce an OH airglow spectrum, or by using a combination of noble gas 
arc lamp spectra. The wavelength calibration image is spatially and spectrally rectified onto a 
rectilinear grid using predetermined maps of the distortion in each of these directions. A cubic 
function of fixed cubic coefficient and zero quadratic coefficient (determined empirically to fit the 
distortion well) is fit to the wavelength reference, providing a mapping from pixel coordinates to 
linear wavelength coordinates. This mapping is combined with the predetermined distortion maps, 
and then applied in a single interpolation to the data frames. Once rectified, the pairs are combined 
into a single spectrum, using Gaussian centroids for the spatial registration, or known pixel offsets 
for galaxies having diffuse or faint continua. This process removes any residual sky subtraction 
signal, as one effectively subtracts the simultaneous sky at the alternate slit position in the process. 
The noise at each pixel in the two-dimensional spectrum is computed along with the spectral data 
by quadrature combination of read noise and Poisson shot noise based on signal level for each pixel. 
Pixel replacement, atmospheric calibrator division, and interpolation in the data frames are also 
properly treated in the noise frames throughout the process. The same extractions applied to the 
spectral data are performed on the noise image, producing a one-dimensional noise array along 
with the extracted spectrum. 

The seeing at the time of observation was assessed by taking short guided exposures — typically 
10 s — of nearby field stars through the wide-open slit of the spectrograph, with the spectrograph 
operating in imaging mode. In this way, the point spread function (PSF) of the telescope-plus- 
atmosphere during the spectral observations could be estimated. These PSF exposures were typi- 
cally taken 3-4 times during an observation, usually at the beginnings and ends of the 30-40 minute 
intervals spent at each grating setting. The full-width at half-maximum (FWHM) along the slit 
direction was measured for each PSF image, and the nuclear extractions of the spectra were based 
on this measurement. In each case, the PSF FWHM was rounded up to an integral number of pix- 
els, and extractions matching these integer widths were performed. The centers of these apertures 
were not constrained, however, such that fractional pixels were summed in the production of the 
nuclear spectra. The widths of the extractions, in arcseconds, are listed in Table 2. 
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Table 2. Observational Parameters 



Galaxy 


Grating 


Date 


Integration 


Slit 


Atmos. 


Obj./Cal. 


Extraction 




Setting 




(sec) 


P.A. 


Calib. 


Airmass 




IRAS 00262+4251 


Paa 


3 Jan 1996 


1800 


0° 


G1.5V 


1.02/1.05 


1.5 




[Si VI] 




1800 






1.11/1.04 


1.7 




H 2 +Br7 




1800 






1.05/1.05 


1.5 


IRAS 01521+5224 


Paa 


2 Jan 1996 


1800 


26° 


G5Vb 


1.06/1.09 


2.0 




[Si VI] 




1800 






1.06/1.08 


1.8 




H 2 +Br7 




1800 






1.07/1.08 


1.7 


IRAS 04232+1436 


Paa 


26 Nov 1996 


1800 


5° 


GOV 


1.07/1.04 


0.7 




[Si VI] 




1800 






1.07/1.06 


0.8 




H 2 +Br7 




1800 






1.06/1.06 


0.7 


IRAS 05246+0103 


Paa 


5 Dec 1996 


1200 


109° 


G1V 


1.19/1.18 


0.8 




[Si VI] 




1800 






1.21/1.21 


0.8 




H 2 +Br7 




1800 






1.18/1.17 


0.7 


IRAS 08030+5243 


Paa 


31 Dec 1995 


1800 


151° 


G5V 


1.07/1.07 


1.5 




[Si VI] 




1800 






1.07/1.06 


1.3 




H 2 +Br7 




1800 






1.06/1.07 


1.3 


IRAS 08311-2459 


Paa 


8 Apr 1996 


1800 


0° 


G3V 


1.91/1.99 


1.7 




[Si VI] 




1800 






1.92/1.82 


1.7 




H 2 +Br7 




1200 






1.94/1.93 


1.0 


IRAS 08344+5105 


Paa 


3 Jan 1996 


1800 


0° 


GOV 


1.09/1.10 


1.7 




[Si VI] 




1800 






1.05/1.09 


1.5 




H 2 +Br7 




1800 






1.07/1.09 


1.5 


IRAS 08572+3915 


Paa 


25 Nov 1996 


2400 


150° 


GOV 


1.01/1.01 


0.8 




[Si VI] 




1200 






1.03/1.02 


0.7 




H 2 +Br 7 




1800 






1.01/1.01 


0.7 


IRAS 09061-1248 


Paa 


28 Nov 1996 


1800 


21° 


GOV 


1.49/1.51 


1.5 




[Si VI] 




1800 






1.47/1.43 


2.0 




H 2 +Br7 




1800 






1.45/1.46 


1.8 


IRAS 09111-1007 


Paa 


11 Feb 1997 


1800 


40° 


G2V 


1.39/1.45 


1.7 




[Si VI] 




1200 






1.42/1.39 


1.5 




H 2 +Br7 




1800 






1.39/1.41 


1.7 


IRAS 09583+4714 


Paa 


3 Jan 1996 


1800 


56? 5 


GOV 


1.03/1.03 


1.5 




[Si VI] 




1800 


56? 5 




1.03/1.04 


1.2 




H 2 +Br7 


31 Dec 1995 


1800 


125° 




1.04/1.03 


0.8 


IRAS 10035+4852 


Paa 


22 Feb 1997 


1200 


33° 


GOV 


1.04/1.05 


1.0 




[Si VI] 


23 May 1997 


1800 




G1V 


1.21/1.15 


0.8 




H 2 +Br7 




1800 




F9V 


1.09/1.06 


0.7 


IRAS 10190+1322 


Paa 


22 May 1997 


1800 


65° 


GOV 


1.15/1.16 


1.0 




[Si VI] 




1800 




G1V 


1.26/1.23 


1.2 




H 2 +Br7 




1800 




GOV 


1.44/1.42 


1.0 


IRAS 10494+4424 


Paa 


22 Feb 1997 


1800 


148° 


G2V 


1.03/1.02 


1.0 




[Si VI] 




1200 




GOVs 


1.15/1.10 


0.8 




H 2 +Br7 




1800 






1.09/1.10 


0.7 


IRAS 11095-0238 


Paa 


6 May 1996 


1800 


38° 


G1V 


1.24/1.26 


1.0 




[Si VI] 




1800 






1.29/1.31 


1.0 




H 2 +Br 7 




1800 






1.25/1.29 


0.8 


IRAS 12112+0305 


Paa 


5 May 1996 


1800 


34° 




1.09/ ■ ■ • 


1.0 




[Si VI] 




1800 






1.05/ ■ ■ • 


0.8 




H 2 +Br7 




1800 






1.07/ ■■• 


1.0 
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Table 2 — Continued 



Galaxy 


Grating 


Date 


Integration 


Slit 


Atmos. 


Obj./Cal. 


Extraction 




Setting 




(sec) 


P.A. 


Calib. 


Airmass 


(") 


TR A G 1/1Q/10 1/1/17 
ITrtAO 14o4o — 144 / 


Poo, 


o iviay iyyo 


1 onn 

loUU 


QO° 
OZ 


p 1 v 


1 ^1 /l /IV 
1.01/ 1.4 ( 


i n 

l.U 




[Ol V1J 




1 onn 
loUU 




P/1V 


1 70 / I fiQ 
1. ( Z/ l.Do 


i n 

l.U 




H2+Br7 




1800 






l.Oy/ l.DZ 


1.0 


TP A G 1/1Qc;Q mcd 
lltAo 14oOZ— iy04 


Pact 


o May iyyo 


i onn 
loUU 


no 
U 




1 7n /l 7Q 
1. (()/ 1. IO 


1 o 
l.Z 




rc; VTl 
[Ol V1J 




i onn 

loUU 




KjtU V 


1 71/1 70 

i. < 1/ 1. ( y 


n o 
U.o 




IT 1 ~D „_ 

Il2+-t>r7 




1800 




KsrO V 


1 CO /l 70 

i.oo/ 1. r z 


0.8 


111 /AO 14oy4-t-OOOZ 




o iviay iyyo 


i son 

loUU 


OD 


P 1 V 
ul V 


1.04/ l.OU 


U.o 




rc; VTl 
[Ol V1J 


o iviay iyyo 


i onn 

loUU 






1 k i /i en 
1.01/ l.OU 


n o 
u.o 




H2+Br7 




1800 






1 oc /i on 

l.oo/l.oz 


0.7 


lltAo 10z40+lUly 


Paa 


zz iviay iyy ( 


i onn 
loUU 


1 1 Q° 
llo 


\jo V 


1 11/1 1 o 
1.11/ 1.1Z 


1 o 
l.Z 




Tg; VTl 
[ol V1J 




1 onn 
loUU 






i nn /i 1 1 
i.uy/ 1.11 


i n 
l.U 




TJ 1 ~D „_ 




1800 






1.10/1.10 


1.0 


TP A G 1 ^O^njQfinG 
lltAo lOZOU+oDUy 


"Port, 


zi iviay iyy ( 


1 onn 
loUU 


loO 


POV 


1 nQ /i m 
l.Uo/ 1.U1 


1 n 

l.U 




rc; VTl 
[Ol V1J 




O/i nn 

Z4UU 






i m /i no 
1 .Ul / 1 .uz 


1 n 
l.U 




H2+Br7 




i onn 
loUU 






i nn /i m 
l.UU/ 1.U1 


n 7 
U. / 


TP A G 1 K/lfiO n/i t;n 
lltAo 104DZ — U40U 


Paa 


i Apr iyyo 


1 onn 
loUU 


cn° 

ou 


LjZ.O V 


1 QO /l Q C 

L.oZf l.ou 


1 n 
l.U 




TG; VTl 
[Ol V1J 


o Apr iyyo 


i onn 
loUU 






i /ii /i q^; 
1.41/ l.ou 


n o 
U.o 




H2+Br7 


*7 A -1 n C\ C 

7 Apr 1996 


1800 






1.39/1.37 


1.0 


TP A G 1 «A 07_l_ f^/1 /I 7 
IrlAo lD4o/H-044( 


"Port, 

.r act 


zi iviay iyy i 


1 onn 

loUU 


71 
i 1 


P tvv 

LtO V 


1 no /i 1 1 
l.Uo/ 1.11 


n 7 
u. / 




Tg; vtI 

[Ol V1J 




i onn 

loUU 






1 11/1 1 o 
1.11/1 . iz 


n o 
U.o 




H2+Br7 




1800 






i iv/i on 
1.17 / l.ZU 


n 7 
U. ( 


TP Ad 1 7nOO_l_^01 7 
mAo l/UZo-t-Ool( 


r^aa 


zz iviay iyy < 


1 onn 

loUU 


yo 


LtZ V 


1 1/1/1 1 q 

1 . 14/ 1 . lo 


n o 
U.o 




rc; VTl 
[Ol V1J 




i onn 
loUU 






1 10/1 1/1 

1. lo/ 1.14 


i n 
l.U 




H2+Br7 




1800 




LrU Va 


1 O /I / 1 oc 

1.Z4/ l.zu 


n 7 
U. ( 


TP A G 1 fi/17n_l_QOQQ 
mAo lo4(U-roZoo 


Pact 


Q1 Tnl 1 0Qf; 

oi jui iyyo 


1 onn 

loUU 


P.7° 
D / 


p sv 


1 nc: /i nn 
1.U0/ l.UU 


n o 
U.o 




Tg; vtI 

[Ol V1J 




1 onn 

loUU 






i nn / 1 nn 
l.UU/ l.UU 


1 n 

l.U 




H2+Br7 




1800 






i m /i nn 
1.1)1/ l.UU 


1.0 


T"DAQ 1 1 HO A A 

lKAo ly4oo+Uy44 


Paa 


O A n 1 finC 

z Aug iyyo 


1800 


121° 


r OV+CjjUV 


i nn / 1 nn 

i.uy/i.uy 


0.7 




Tg; vtI 
[Ol V1J 




i onn 
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The most significant stellar absorption feature present in the 2 fj,m window for the atmo- 
spheric calibrator stars is the Br7 line at 2.1655 /mi. While a template stellar spectrum may be 
adequate for removal of this line, we instead interpolated over the line in the two-dimensional G 
star calibration/flat-field spectrum before dividing the galaxy data by the calibrator data. Multi- 
plying the corrected galaxy spectrum by a Planck blackbody function with a temperature matching 
that of the calibrator rectifies the spectrum to a proper representation in flux density. To remove 
the effects of the weaker absorption lines in the atmospheric calibrator, the galaxy spectrum was 
multiplied by a normalized solar spectrum, smoothed and re-binned to the same resolution. For 
the long wavelength grating setting covering the H2 1-0 S(l) and Br7 lines, the CO band structure 
becomes significant when using late-type G stars for the atmospheric calibration. When the atmo- 
spheric calibrator was of type G4-G8, a linear combination of G3 V and K0 V template spectra 
from Kleinmann & Hall (1986) was used instead of the solar spectrum. Again, this template was 
smoothed and re-binned to match the individual galaxy spectra. 



4. Results 

The chief results from this spectroscopic survey — that AGN appear to be rare among ULIRGs — 
has been reported previously in Murphy et al. (1999). In the present paper we focus attention on 
the range of observed properties rather than on the collected sample as a whole. 

Figure 1 presents spectral extractions for the primary nuclei in the 33 sample galaxies, with 
aperture sizes indicated in Table 2. Each plot is fixed in the rest frame, with the observed wavelength 
scale indicated at bottom. Annotations mark the locations of every common emission line, though 
in any individual spectrum only the lines with reasonable detections are labeled. The location of the 
tick marks is based on the observed redshift of the Pact line, as is the rest wavelength scale at top. 
Each spectrum is composed of three independent spectral datasets. For each of the three spectral 
ranges, a best linear fit to the line-free continuum is assessed for the purpose of determining an 
overall average slope for the spectrum. The individual spectra are then scaled to lie on a common 
linear "backbone" defined by the average slope, and combined in a noise-weighted manner across 
the overlap region around 1.92 /im. The vertical scale is normalized such that the continuum is at 
one unit of flux density at 2.155 //m, corresponding to the middle of the K s bandpass. 



Table 2 — Continued 
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Figure 1: continued 



Figure 2 contains spectra of the secondary nuclei in the 14 sample galaxies with a nearby 
galactic companion, employing the same wavelength scale used for Figure 1. Annotations and 
scalings are the same as for Figure 1, though some of the weaker lines are not marked by default 
in this set. The tick marks denoting line positions are based on the redshift of the Paa line in 
the primary nucleus, so that velocity displacements between the two nuclei are more easily seen. 
The numerical flux density scale for each is identical to that for the corresponding spectrum in 
Figure 1, allowing direct comparison of the relative levels for each. The relative scaling of the three 
individual spectra determined for Figure 1 are used in Figure 2 as well. 

The choice of primary versus secondary nucleus is based solely on the continuum level, with 
the brighter nucleus designated as the primary. The only exception is IRAS 22491—1808, for 
which the primary nucleus is identified with the dominant line-emitting, yet weaker continuum 
nucleus to the east. The western nucleus is 30% brighter, but lacks the strong H2 emission that is 
associated with every other primary ULIRG nucleus. Comparison of Figures 1 and 2 indicate that 
the primary nuclei, besides being the brightest continuum sources, also exhibit apparently stronger 
line emission. Yet despite the weaker line strengths, the secondary nuclei often do show appreciable 
line emission activity, appearing relatively similar to the primary nuclei. If the observed Paa line 
strength is proportional to star formation rates — and therefore total luminosity — then even though 
the dominant luminosity can be associated with one nucleus, the lesser nucleus still contributes a 
substantial fraction of the total luminosity. It must be realized, however, that the high levels of 
extinction present in ULIRGs could easily lead one to misjudge the identity of the primary source 
of infrared luminosity based simply on apparent line strengths. 

Table 3 presents the line measurements from the primary nuclear data shown in Figure 1 and 
Table 4 presents the same data for the secondary nuclear spectra of Figure 2. For each galaxy, 
the continuum slope and measured Paa FWHM are given, as well as equivalent width, FWHM, 



-19- 



rest wavelength (urn) 

2.00 2.10 



2.00 



2.1 



2.20 



2.30 



2.20 




2.40 



1 .80 


1 .90 


2.00 2.10 2.20 


3 1.2 






RAS 08572+3915SE 




1 0.8 




8 






CD 

x 0.6 

u 




D 
Q_ 






^ 0.4 

a> 










'1 0.2 
^ 0.0 








n 



1 .80 

2.0 

1 .5 - 

1 .0 - 

0.5 
0.0 



2.00 
1 .90 



2.10 
2.00 



2.20 



2.30 



2.10 



2.20 



>> 

en 
c 

CD 



CD 
> 



_o 

CD 




RAS 09061 -1 248N 





2.00 



2.10 2.20 
observed wavelength (/xm) 



2.30 



-20- 



1 .80 
3 2.5 — 

2.0 - 



to 
c 

a; 1 

T3 



CD 
> 



CD 



5 b- 
.0 
.5 



0.0 



1 .90 



rest wavelength (urn) 

2.00 2.10 



2.20 



I 



RAS 09583+4714SW 





1 .80 



CO 
£_ 
C 

~o 

X 

"1 



> 



^3 

CJ 




2.00 



2.1 



2.20 



2.30 




_>N 

'to 

c 
c 

T3 



C 

> 



5 f- 

4 §- 

3 |- 

2 
1 

rL 




RAS 121 12 + 0305NE 



2.00 



2.10 2.20 
observed wavelength (/xm) 



2.30 



- 21 - 



1 .80 
4 



c 

CD 

T3 



CD , r- 

> I ~ 

CD 



1 .90 



rest wavelength (um) 

2.00 2.10 

RAS 14348-1447NE 



Paa 



CD 



no 



2.20 




en 
c 

CD 
"D 

X 
Z3 



CD 
> 



CD 



CO 
C 
CD 
T3 



CD 
> 



O 
CD 



'en 

c 

CD 



CD 
> 



CD 




2.00 



2.10 2.20 2.30 

observed wavelength (/xm) 



2.40 



-22- 




Figure 2: continued 
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and velocity centroid of each significant line. In this context, weaker hydrogen recombination 
and molecular hydrogen vibration-rotation lines are excluded from the table as their properties 
follow from those of the stronger lines. The continuum slope, a, is reported as the power-law 
index, where /a(A) = C\ a . The equivalent linear slope, with f\(\) = /a(Ao) [1 + c(A — Ao)], where 
Ao = 2.155 jum, is a = a/\ . The measured line properties reported in Tables 3 and 4 are obtained 
by the procedures detailed in Appendix A. 



Table 3. Measured Line Properties for Primary Nuclei 
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377+38 
-2+29 

5.65+0.25 
350+49 
-59+26 

1.095+0.086 
295+63 
-65+42 

3.15+0.12 
450+35 
-65+35 

1.19+0.11 
193+58 
9+37 

5.03+0.097 
490+22 
-5+31 

1.324+0.067 
307+46 
-19+30 

1.45+0.14 

323+88 
6+35 

0.48+0.029 
297+53 
229+40 

7.13+0.29 
374+36 
-52+33 

1.26+0.078 
250+48 
3+33 

2.04+0.16 
305+68 
11+44 



0.154+0.072 
0+130 
-50+140 



0.34+0.14 
220+290 
-100+190 

0.244+0.083 
120+210 
170+110 



0.143+0.08 
160+280 
210+170 

0.341+0.078 
300+210 
137+67 

0.55+0.1 
400+180 
-15+60 



0.35+0.13 
140+170 
120+130 

0.41+0.16 
200+300 
140+120 



1.248+0.096 
336+65 
7+32 

1.89+0.11 
427+65 
-3+29 

8.56+0.61 
241+43 
-73+35 

0.957+0.071 
0+31 
-19+30 

3.11+0.14 

334+36 
-46+30 

1.88+0.15 

243+48 
35+34 

4.9+0.16 
426+33 
61+33 

1.242+0.088 
247+50 
25+27 

1.928+0.059 
223+19 
58+32 

0.268+0.024 
102+43 
158+32 

5.95+0.39 

265+54 
-134+30 

1.53+0.14 
301+74 
6+39 

2.46+0.17 
319+49 
11+34 



1.12+0.11 

320+91 
22+34 

1.73+0.16 
364+75 
-51+31 

2.5+1 
0+100 
80+130 

1.57+0.15 
185+53 
-51+33 

1.35+0.2 
127+80 
35+40 

2.97+0.26 
166+46 
6+32 

2.14+0.27 
370+110 
127+47 

2.08+0.13 
304+51 
-15+25 

1.429+0.077 
204+34 
100+33 

0.516+0.06 
360+110 
237+44 

3.46+0.89 
230+140 
-195+78 

2.13+0.25 
336+94 
22+45 

3.52+0.25 
249+50 
-101+30 
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Figure 3 presents the two-dimensional spectra around the Paa line for each of the sample 
galaxies. In each, the horizontal axis origin is located at the center of the primary nucleus. The 
apertures corresponding to the extractions displayed in Figures 1 and 2 are indicated at the top 
and bottom of each panel. Annotations to the left and right of each panel denote the sense of 
the slit orientation, the exact position angles for which are listed in Table 2. The two-dimensional 
spectrum for IRAS 09583+4714 has been split into two panels, though the horizontal scale reflects 
the real separation between the two nuclear components. 

4.1. Hydrogen Recombination & Extinction 

The K band spectra of the galaxies in this sample are dominated by the Paa line, which 
is the strongest line accessible to search for signs of broad emission. Only in two of the sample 
galaxies, IRAS 08311—2459 and IRAS 15462—0450, is any such emission apparent. Both of these 
galaxies are optically classified as Seyfert galaxies, so it seems that obscured AGN in ULIRGs are 
either extremely rare, or suffer a few magnitudes of extinction even at near-infrared wavelengths, 
as was also suggested by Goldader et al. (1995). A few galaxies, such as IRAS 14348—1447 and 
IRAS 14394+5332 have excess emission on the blue side of the Paa line, but in both cases there 
appears to be emission from He I contributing significantly to this flux, and no appreciable emission 
on the red side of the narrow Paa profile. Furthermore, the H2 lines in IRAS 14394+5332 have 
an asymmetric blue shape more typical of wind or outflow phenomena than of velocity-broadening 
due to an AGN. IRAS 20087—0308 has some hint of broad emission at the base of the Paa line, 
though not significantly elevated out of the noise, and again confused by He I. As seen in Figure 3, 
the moderate 600 km s _1 width of the narrow portion of the Paa profile in IRAS 20087—0308 can 
be attributed to a steep rotation gradient on the nucleus. The apparent broad base on the Paa 
line in the southwestern nucleus of IRAS 12112+0305 may be spurious, as this is the one galaxy 
in the sample whose spectra were not calibrated by a G star observation, but rather by a template 
atmospheric spectrum. The excess blue emission happens to be coincident with the significant CO2 
atmospheric absorption at 2.00-2.02 jtxm, so that a slight error in estimated column depth results 
in a large localized continuum offset in the corrected spectrum. 

The two-dimensional spectra of Paa in Figure 3 demonstrate that the atomic recombination 
line emission in ULIRGs is almost always dominated by a nuclear component. Several ULIRGs show 
off-nuclear emission at lower levels, and a couple of sources appear to be primarily characterized 
by diffuse line emission — namely IRAS 04232+1436 and IRAS 20046—0623. When two nuclei are 
present, the secondary nucleus generally shows emission in Paa as well. Exceptions to this are 
IRAS 16487+5447 and IRAS 23327+2913. The extra continuum object seen in the panel with 
IRAS 08030+5243 is a nearby star. 

The fraction of the total line emission that is unresolved can be estimated for each ULIRG by 
computing the flux contribution from a Gaussian spatial profile having the width of the estimated 
seeing and sharing the maximum amplitude of the spatial line emission distribution. Comparing 



Table 3 — Continued 



Galaxy 


Continuum Paa PWHM b 


Rows c 


He I Pact 


H 2 S(3) 


[Si VI] [Fc II] 


H 2 S(l) 


Br7 




Slope" (km s- 1 ) 




1.8689 1.8751 


1.9570 


1.9629 1.9670 


2.1213 


2.1655 


IRAS 19458+0944 


-0.25 269 


















EW 


0.668+0.086 11.093+0.075 


1.46+0.13 




1.236+0.089 


1.74+0.29 






FWHM 


390+180 272.9+5.1 


324+84 




237+50 


150+110 






cz 


-27+58 29982+35 


-83+59 




3+47 


73+53 


IRAS 20046-0623 


-3.99 244 


















EW 


6.3+0.21 


1.16+0.2 


0.45+0.2 


1.63+0.44 


3.67+1 






FWHM 


247+27 


190+100 


0+140 


0+88 


310+170 






cz 


25293+14 


-16+43 


50+140 


6+60 


100+80 


IRAS 20087-0308 


-0.3 595 


















EW 


0.465+0.093 10.22+0.14 


1.611+0.067 




2.11+0.15 


0.64+0.24 






FWHM 


270+150 574+23 


475+49 




570+110 


280+280 






cz 


-92+48 31600+20 


-41+33 




64+30 


36+97 


IRAS 20414-1651 


-1.93 515 


















EW 


0.82+0.16 10.47+0.21 


1.64+0.16 




1.71+0.17 


2.23+0.63 






FWHM 


400+190 488+20 


287+82 




460+120 


170+180 






cz 


53+71 26043+25 


16+49 




35+46 


-137+92 


IRAS 21504-0628 


-1.29 256 


















EW 


l.Uf±0.14 lo.o2±U.14 


U . o * ± U . 1 4 




i.oo4ztu.uyy 


1 V/i _i_n 1 q 
1. * 4±U. lo 






FWHM 


230+100 260.8+5.6 


250±130 




310±43 


144±40 






cz 


-17+52 23352+26 


46 ±44 




43 ±43 


-4±41 


IRAS 22491-1808E 


-4.12 128 


















EW 


13.21+0.38 


3.06±0.23 




b.26±0.3 i 


rf.U4±U.»0 






FWHM 


131 + 18 


0+52 




302+60 


430±350 






cz 


23277+24 


-35+41 




-2+33 


10±77 


IRAS 23327+2913S 


0.08 151 


















EW 


0.43+0.11 6.85+0.12 


0.87+0.26 




1.92+0.36 








FWHM 


170+130 193+9 


110+120 




210+120 








cz 


68+76 32235+18 


-1+77 




219+51 




IRAS 23365+3604 


-1.46 186 


















EW 


9.284+0.083 


1.291+0.079 




1.493+0.057 


1.345±0.094 






FWHM 


191.1+5 


245+45 




219+22 


216±39 






cz 


19309+17 


60+24 




-36+25 


-9±28 


a Slopc expressed a 


s power-law index, a, where fx — 


c\ a . 












b Strictly formal measure of the full-width at half-maximum of the Paa line, deconvolved by the 


instrumental r 


csolution. 






c For each galaxy, 


the following measures are provided for selected lines: equivalent width (EW) is given in 


nanometers; FWHM as 


estimated via li 


nc flux and peak 


amplitude (sec text), 


deconvolved by the instrumental 


resolution; 


first moment ccntroid of the line, 


expressed as 


recessional velocity (cz) for Pact, and vc] 


ocity offsets from 


Pact for all other lines. The first moment measure looses accuracy 


for weak lines depending on the spectral range 


summed, and therefore 


is of limited use 


here. 



to 



Table 4. Measured Line Properties for Secondary Nuclei 



Slope" Pact PWHM b Rows c Ho I Pac* H 2 S(3) [Fe II] H 2 S(l) Br 7 



Galaxy 



IRAS 01521+B224N -4.18 



IRAS 05246+0103W 2.54 



IRAS 08572+3915S 7.02 



IRAS 09061-1248N -7.07 



IRAS 09583+4714SW -6.81 



IRAS 10035+4852SW -10.31 



IRAS 10190+1322W -3.94 



IRAS 12112+0305NE 2.53 



IRAS 14348-1447NE -2.11 



IRAS 15245+1019W -5.05 



IRAS 17028+5817E -1.31 



IRAS 22491-1808W -8.88 



EW 
FWHM 



EW 
FWHM 



EW 
FWHM 



EW 
FWHM 



EW 
FWHM 



EW 
FWHM 



EW 
FWHM 



EW 
FWHM 



EW 
FWHM 



EW 
FWHM 



EW 
FWHM 



EW 
FWHM 



0.46±0.16 

0±140 
-60±120 



0.59±0.23 
50±160 
-60±130 



2.06±0.3 
450±160 
133±64 



0.93±0.38 
0±140 
80±100 



1.12±0.6 
190±230 
23990±170 

7.93±0.31 

72±19 
29103±23 

2.92±0.21 

0±19 
17465±22 

5.52±0.59 

252±76 
22030±30 

1.97±0.16 

246±41 
25696±31 

10.33±0.54 

430±31 
19443±31 

4.51±0.14 

271±26 
23012±21 

21.49±0.52 

271±19 
21875±20 

15.41±0.35 

272±14 
24681±23 

4.36±0.38 

310±70 
22577±31 

9.66±0.76 

97±33 
31938±32 

5.59±0.2 

0±23 
23274±23 



0.84±0.29 

0±170 
-60±140 

0.94±0.19 
320±180 
138±76 

1.03±0.13 
320±110 
-61±54 

0.41±0.092 
40±110 
-70±64 

4.62±0.13 
401±27 
-110±36 

1.98±0.15 
218±63 
-31±35 

2.7±0.19 

406±74 
-72±42 



0.52±0.21 

0±99 
-180±160 



0.31±0.11 

0±150 
370±170 



0.94±0.29 
60±100 
126±77 



2.13±0.39 
253±95 
81±61 



1.88±0.14 
292±47 
26±42 

0.64±0.13 
63±89 
47±56 

4.05±0.2 
366±49 
-47±30 

2.5±0.25 
274±76 
-66±36 

2.19±0.24 

78±54 
-100±41 



1.36±0.17 
212±67 
-17±48 

0.72±0.18 
150±120 

-28±68 

3.67±0.39 
278±78 
-12±36 

1.64±0.5 
0±120 
20±100 



Note. — 



Sec notes to Table 3 for details. 
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Fig. 3. — Two dimensional spectra around the Paa line for each galaxy. The horizontal scale is in 
arcseconds, with the primary nucleus at the origin. The spectral direction spans about 11 resolution 
elements, or about 3300 km s -1 . Extraction apertures are indicated on the top and bottom of each 
frame. 
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Figure 3: continued 
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Figure 3: continued 
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this flux with the measured flux provides an upper limit to the fraction of light contained within 
an unresolved core. Most of the primary ULIRG nuclei are consistent with being unresolved, with 
a median seeing of l'.'O and a median physical resolution of 1.4 kpc. The scatter in the estimated 
fraction of unresolved flux is at the ±20% level, mostly due to uncertainty in the seeing estimate. 
Only 8 of the 33 ULIRGs have estimated unresolved contributions lower than 80%, such that the 
majority of ULIRGs contain most — if not all — of their Paa emission within the central 1 kpc radius. 
Only three ULIRGs in the sample (IRAS 04232+1436, IRAS 18470+3233, and IRAS 20046-0623) 
have less than 50% of their Paa flux in an unresolved core. 

The sample was chosen to allow coverage of both the Paa and Br7 emission lines so that a 
comparison between these lines could be made in order to assess extinction to the line emitting 
source. Table 5 lists the inferred visual extinction values for the sample galaxies. These estimates 
assume an intrinsic ratio for Paa/Br7 of 12.1 (Osterbrock 1989, Case B with n e = 10 4 cm -3 and 
T = 10000 K), and Ap aa = 0A45Ay and Asr-y = 0.116^4y (interpolated from the extinction law 
of Rieke & Lebofsky 1985). 

Though there exists very little leverage in wavelength between Paa and Br7, some ULIRGs are 
apparently obscured enough to show definite reddening across this short range. Typical la errors 
on the extinction estimates are around 5 magnitudes, though the actual errors could be larger if 
disparate aperture sizes were used in the spectral extractions. Indeed, each of the four ULIRGs 
with extinctions more negative than la have aperture differences between the Paa and Br7 spectra 
of at least 50%. Almost half of the ULIRGs have estimated extinctions greater than the la limit, 
and 11 have at least 2a extinction estimates — typically corresponding to Ay > 10 mag, with many 
in the range of Ay ~ 25 mag. 

Figure 4 presents the extinction measurements in histogram form. Different shading represents 
varying levels of confidence in the measurement. Higher extinction estimates are allowed greater 
variance in classifying their reliability, provided they are significantly different from zero extinction. 
Counting only the 22 ULIRGs with reliable extinction measures, the median visual extinction of 
~ 8 mag corresponds to an extinction at Paa of approximately one magnitude. This is a lower limit 
to the nuclear extinction because we can only estimate extinction to the gas we actually detect at 2 
/im, which very likely is peripheral to the nuclear region. Nonetheless, one magnitude of extinction 
at Paa would not readily hide the presence of a central AGN. The very large visual extinctions 
of 25 mag seen in some ULIRGs translate to almost 4 mag of extinction at Paa, which almost 
certainly would hide an AGN from view with the currently achieved signal-to-noise ratio in the 
ULIRG spectra. 

Because an active galactic nucleus is believed to indicate the presence of a QSO (quasi-stellar 
object) in the galaxy nucleus, one may quantify the amount of obscuration necessary to hide a 
central AGN by estimating the extinction required to hide the line flux expected from the high- 
velocity wings of Paa in a typical QSO of the appropriate bolometric luminosity. For this purpose, 
we use the median QSO spectrum presented in Murphy et al. (1999) to represent the appearance of 
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Table 5. Derived Properties 



Galaxy Extinction a H2 Line Parameters 







T ro t 


T-vib 


H2 Excitation 




(Av) 


(K) 


(K) 


Mechanism 


IRAS 00262+4251 


+ 66 


1300 ± 75 

1 WWW — 1 — 1 KM 


4000 ± 1000 

I WWW — 1_ 1 W V ' w 


mixture 


IRAS 01521+5224 


3 + 14 


1400 ± 300 

1 I WW — 1 — • ) WW 


5000 


thermal* 


IRAS 04232+1436 


2.3 ± 7.2 


2100 ± 225 




thermal 


IRAS 05246+0103 


1 + 15 


2500 ± 375 

>\ )\ ) — 1 — • ) t w 




thermal 


IRAS 08030+5243 


0.9 ± 4.6 


2100 ± 450 

ilUU — 1 — 1 w w 




thermal 


IRAS 08311—2459 


— 14.5 ± 4.3 d 


4300 ± 500 

1 W W W — 1 J WW 




thermal 


IRAS 08344+5105 

1 1 V - i i j V 7 7' * 7 1 1 • y 1 WW 


—5.0 ± 5.8 


1800 ± 200 

iuuu — i— — v / w 


4500 ± 2000 

1 w w w — 1_ ^- w w w 


mixture 


IRAS 08572+3915 


12.1 ± 3.6 


2500 ± 300 

£j%j\J\j — 1 WWW 


5000 ± 1500 

W W W W — 1_ 1 J WW 


mixture 


TRAS 09061—1248 

11 lii kJ \J JUU -1 i w 71 W 


3.9 ± 6.1 


2200 + 250 

Z- Z- W W 1 Z- t. ) \ I 


3500 


IjllWl 111 CXI 


TRAS 091 1 1 —1007 


8.7 ± 4.0 


2300 + 200 

W W W i Z_ V / V / 


2700 


fViprin a 1 * 

L11C1 111 ell 


TRAS 09583+4714 


—6.2 ± 5.1 d 


2000 + 250 




fViprinfll 

Liir^j 111 *xi 


IRAS 10035+4852 


—9.9 ± 5.3 d 


2900 ± 300 

£-t\J\J\J 1 WWW 




thermal 


IRAS 10190+1322 


4.4 ± 3.9 


2200 ± 125 




thermal 


IRAS 10494+4424 


— 10.8 ± 3.4 d 


2400 ± 250 


3500 


thermal* 


TRAS 1 1095—0238 


9.5 ± 16 


1 700 + 1 50 

J- ( WW 1 J- WW 


3800 + 1 300 

<J W W W 1 ± w w w 


TTii^rtiifp 

1 JlljV Lr HI \Z 


IRAS 12112+0305 


19.5 ± 3.7 


2500 ± 300 

W W W — 1 WWW 


5000 


t hernial* 


IRAS 14348—1447 


5.7 ± 5.7 


2000 ± 125 

£j\J\s\J — 1 1ZU 




thermal 


TR AS 14359— 1 954 

J_J.ljiT.lJ ± t: » J W J. £/W7t 


2.8 ± 3.2 


1 600 + 1 50 

1UUU _l_ lw 




+ n r*T""m 'A 1 

I'llLl lllcll 


TRAS 14394+5332 


12.0 ± 4.8 


2200 + 1 50 

iiZUU 1 1UU 


2800 


fViprm fll * 

LlltM 111 ell 


IRAS 15245+1019 


7.3 + 2.3 


4900 ± 800 




thermal 


IRAS 15250+3609 


4.6 + 2. l d 


2000 ± 150 


3100 ±600 


mixture 


IRAS 15462-0450 


8.2 + 4.4 


> 4000 




thermal 


IRAS 16487+5447 


23.3 + 9.9 


1700 ± 150 




thermal 


IRAS 17028+5817 


27.2 + 4.4 


2200 ± 200 




thermal 


IRAS 18470+3233 


0.5 + 2.7 


1600 ± 100 




thermal 


IRAS 19458+0944 


24.3 + 6.3 


2300 ± 100 




thermal 


IRAS 20046-0623 


50 + 11 


1900 ± 300 




thermal 


IRAS 20087-0308 


-10 + 15 


2200 ± 350 




thermal 


IRAS 20414-1651 


26 + 11 


2400 ± 225 




thermal 


IRAS 21504-0628 


3.4 + 2.8 


2600 ± 300 




thermal 


IRAS 22491-1808 


16 + 11 


2100 ± 250 




thermal 


IRAS 23327+2913 


2 + 73 


1600 ± 200 




thermal 


IRAS 23365+3604 


14.5 + 2.6 


1900 ± 150 


4200 ± 900 


mixture 



a Based on Paa/Br7 line ratio. 

b Given only for definite thermal/fluorescent mixtures (with errors) and for likely 
mixtures (when both of the observed 2-1 lines are above the rotation temperature 
prediction). The vibration temperature assumes an ortho-to-para ratio of 2:1 in 
the 2-1 transition states. 

c Classified as consistent with purely thermal, definite mixture of thermal and 
fluorescence, or possible mixture of the two (with asterisk). 

d Apertures for the Paa and Br7 spectral extractions differ by at least 50%. 
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an unobscured quasar. Identifying the region around 1.881 /jm (« 1000 km s _1 redward of the Paa 
line center) as being well away from the narrow Paa emission, we measure the median QSO flux 
density to be 22% higher than the QSO continuum level. It is found that the broad Paa emission 
from an exposed QSO contributing half the bolometric luminosity of a ULIRG would be comparable 
at 1.881 /xm to the typical ULIRG continuum emission (neglecting the contribution from the QSO 
continuum). Therefore, hiding the presence of the broad emission requires an attenuation factor 
that is comparable to the signal-to-noise ratio of the ULIRG continuum. In the case of the median 
ULIRG spectrum presented in Murphy et al., this translates to an extinction at Paa of 4.4 mag, and 
therefore a visual extinction, Ay ~ 30 mag. The spectra of the individual ULIRGs presented here 
typically have continuum signal-to-noise ratios around 20, so that Ap aa = 3.1 mag, or Ay = 20- 
25 mag of extinction is required to hide a QSO at half the ULIRG luminosity. Details on this 
calculation appear in Appendix B. 

The estimates of extinction made in this way are conservative in two regards. First, the 
aperture photometry of ULIRGs used to establish typical continuum flux densities are based on 
5" apertures, such that extra flux is included from diffuse emission away from the nucleus. In 
smaller apertures corresponding to the spectral extractions employed in this paper, the ULIRG 
continuum is diminished, while the spatially unresolved QSO component would be hardly affected, 
thus increasing the QSO-to-ULIRG contrast considerably, making the QSO more difficult to conceal. 
The second conservative aspect is that the criterion of forcing the broad Paa from the QSO to lie 
at the la level on the ULIRG continuum would not hide the broad emission very effectively, the 
reason being that the 1<7 level is computed for a single pixel (one-quarter of a resolution element). 
When summed over a few resolution elements (i.e., up to the narrow shoulder of the Paa line), the 
broad signal would stand out at more than the 3cr level. 

The result that, given the current level of sensitivity, at least 25 magnitudes of visual extinction 
are required to obscure a QSO in a typical ULIRG, together with the paucity of detections of low- 
level broad line emission — especially evident in the median spectrum of Murphy et al. (1999) — leads 
us to conclude that AGN are either rare in ULIRGs, or very deeply buried. There are five ULIRGs 
in the present sample with measured extinctions based on the Paa/Br7 ratio that are high enough 
to hide AGN without requiring additional attenuation. We are unable to determine the likelihood 
that the rest of the ULIRGs have the requisite nuclear extinctions to hide the presence of AGN. 

The technique of comparing ULIRG spectra to a typical QSO spectrum can also be used to 
estimate the extinction to the broad- line region seen in IRAS 15462—0405. Based on the flux 
density at 1.881 //m, the Paa emission contributes an additional 9% to the continuum level, some 
of which is likely continuum emission from the QSO itself. Correcting for this, one finds an excess 
over the assumed ULIRG continuum (60% of the total continuum) measuring 15%, suggesting an 
extinction Ap aa = 2.0 ± 0.8 mag, or Ay = 14 ± 6 mag. 
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4.2. Molecular Hydrogen Emission 

The most abundant species of line emission in the ULIRG spectra arise from the vibration- 
rotation spectrum of molecular hydrogen. The molecular hydrogen emission shows up strongly in 
all of the primary nuclei of this sample. Three or four lines are typically present in each spectrum, 
most notably the odd rotation state 1-0 transitions. The H2 line strengths vary by more than one 
order-of-magnitude relative to Paa. Figure 5 shows the range of observed line ratios, comparing the 
H2 1-0 S(l) line to Pao. A typical ULIRG has a H2 1-0 S(l)/Paa ratio around 0.15, varying from 
~ 0.05 to 0.40. One galaxy, IRAS 00262+4251 exhibits H2 emission even stronger than the Paa 
line. The bins in Figure 5 are shaded to reflect the classification of the H2 excitation mechanisms, 
as will be discussed below. 

The apparent weakness of the even rotation transitions is an indication that the emission is 
thermally excited, as the thermal condition produces three times as many odd (ortho) states as 
even (para) states. This so-called ortho-to-para ratio is more typically in the range of 1.0-1.8 for 
fluorescent excitation conditions (Black &; van Dishoeck 1987). In addition, the lack of strong 
emission from the 2-1 vibration transition indicates a thermal origin. 

Population diagrams were constructed for each of the primary nucleus extractions in the sam- 
ple, and each characterized in terms of rotation temperature, vibration temperature (where appro- 
priate) and classified in terms of thermal or mixed thermal/fluorescent excitation. The construction 
of population diagrams is discussed in Appendix C. Examples of these diagrams are shown in Fig- 
ure 6. The even rotation state transitions are plotted as two points representing ortho-to-para 
ratios of 3:1 and 1:1. The 3:1 points, represented by boxes in Figure 6, are almost always preferred 
over the 1:1 points (asterisks) appearing below them. The great majority of ULIRGs in this sample 
have no detections in the 2-1 transition lines, but with error bars consistent with purely thermal 
excitation. Solid lines in Figure 6 represent rotation temperature fits, and dashed lines, where 
present, indicate the vibration temperature separating 1-0 and 2-1 points. A few ULIRGs appear 
to be definite mixtures of thermal and fluorescent excitation, with the population diagrams show- 
ing the 2-1 data points well above the rotation temperature line passing through the 1-0 points. 
A number of other ULIRGs have somewhat elevated 2-1 measures, with both the S(2) and S(4) 
lines above the 1-0 rotation temperature line, but with error bars that are still consistent with 
pure thermal emission. The excitation characterizations are summarized in Table 5. For the mixed 
cases, a vibration temperature estimate is given, assuming that the high energy states exist with 
a lower ortho-to-para ratio, likely around 2:1. This has been observed in other mixed excitation 
cases (e.g., Tanaka et al. 1989), and is explained by the fact that the enhanced 2-1 transition 
emission arises from molecules that are primarily excited by fluorescent processes, and thus reflect 
an ortho-to-para ratio closer to that expected for pure fluorescence. 

One point that should be mentioned in the context of the line population diagrams is the 
possible blend of a 2.0024 jim [Fe II] line with the H2 2-1 S(4) line at 2.0035 /jm. The expected 
flux ratio between the [Fe II] line at 1.9670 fim and that at 2.0024 /im is 2-3:1. With the 1.9670 
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Fig. 4. — Histogram of estimated visual extinctions to the Paa line emission for the 33 primary 
ULIRG nuclei. Different shadings represent varying levels of reliability, with most of the negative 
estimates owing to mismatched seeing conditions and therefore significant extraction aperture dif- 
ferences. The mean and median points are marked for the reliable set of measurements. One point, 
corresponding to Ay ~ 50, is beyond the plot boundary. 
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Fig. 5. — Ratio of the H2 1-0 S(l) line strength to that of Paa for the 33 primary ULIRG nuclei. 
A linear scale appears at top. Typical la errors are comparable to the bin width, though perhaps 
a little smaller. The mean is performed on the linear scale, rather than on the logarithmic scale. 
Shadings indicating excitation mechanisms are included for reference. 
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Fig. 6. — Molecular hydrogen population diagrams for six of the ULIRGs, showing the range of 
properties seen in the H2 excitation. The top two show examples of mixed thermal and fluores- 
cent excitation, while the bottom four are consistent with a single temperature characterization. 
IRAS 10494+4244 and IRAS 14394+5332 may have some slight component of fluorescent excita- 
tion, judging by the elevated 2-1 lines. For even (para) rotation states, two points are plotted, 
corresponding to ortho-to-para ratios of 3:1 and 1:1, with boxes corresponding to the 3:1 case, and 
asterisks for 1:1. The population diagram for IRAS 18470+3233 is fairly typical of most of the 
ULIRGs in this sample in terms of excitation information and signal-to-noise ratio. Line fits are 
accompanied by temperatures corresponding to the slopes. 
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/j,m line appearing to be present in a number of ULIRGs, it is to be expected that the H2 2-1 S(4) 
line will suffer some contamination by the [Fe II] line. Indeed, the population diagrams often show 
an anomalously high 2-1 S(4) line flux, with an incompatibly low 2-1 S(2) flux. There is a good 
correlation between such cases and galaxies for which the 1.9670 /im [Fe II] line is believed to be 
present. 

Figure 7 illustrates the distribution of H2 rotation temperatures, indicating a strong clustering 
around 2200 K. Two of the three ULIRGs with estimated rotation temperatures in excess of 4000 
K are classified as having AGN. It is possible in these cases that a significant fraction of the H2 
emission arises from molecular gas excited by X-ray radiation from the central source. Draine & 
Woods (1990) found this mechanism to be the most viable source of H2 excitation in NGC 6240, 
with the X-rays in this instance most probably arising from frequent supernovae or from merger- 
induced shock activity. Table 8 of Appendix C provides expected H2 line ratios in the thermal 
regime across the 1.8-2.4 /im wavelength range. 

Referring to Figure 5, there does not appear to be an overwhelming trend between H2 line 
strength — relative to Paa — and excitation mechanism, except for a possible slight tendency for 
mixtures to be seen among the brighter lines. This could simply be a consequence of higher signal- 
to-noise levels for the weak 2-1 lines in those galaxies with bright H2 lines. More mixed states may 
be revealed with greater sensitivity, most likely in those cases identified here as possible mixtures. 

4.3. Other Emission Line Strengths 

The He I line 1 is detected in the majority (27/33) of the primary nuclear spectra, and in a 
number of the secondary spectra. The noise-weighted mean of detected He I/Paa is 0.058, with a 
relatively small dispersion of ~ 0.03 about this value. The median He I/Paa ratio is 0.063 among 
detected lines, and 0.056 including non-detections as zero-valued points. The full range of detected 
ratios goes from 0.02 (IRAS 15245+3609) to about 0.11 (IRAS 14348-1447; IRAS 14394+5332). 
Lower resolution near-infrared spectra could easily mistake this line for a blue wing on the Paa 
line. 

The [Fe II] line at 1.9670 /im appears in several spectra, most notably the primary nuclei 
of IRAS 00262+4251, IRAS 09111-1007, IRAS 15245+1019, and IRAS 15250+3609, and the 
secondary nucleus of IRAS 10190+1322. These stronger [Fe II] emitters have [Fe II] /Paa ratios up 
to 0.05, though IRAS 00262+4251, with its weak Paa, has a ratio around 0.12. There are many 
possible line detections near the noise limit that have typical [Fe II]/Paa ratios around 0.01-0.02. 



lr The He I emission line just blueward of Paa is a combination of two lines at 1.86860 and 1.86969 (im, both of the 
same Is4/-ls3d transition, but the former a triplet state and the latter a singlet. As such, the shorter wavelength 
line is three times stronger than the longer wavelength line, with a centroid wavelength of 1.8689 fim. The splitting 
between lines is less than half a resolution element for the present spectra, so that the line combination is treated as 
one. 



-39- 



Lines weaker than 1% of Paa are generally too weak to be detected with confidence. The [Fe II] 
line reported for IRAS 08311—2459 is possibly confused with the wing of [Si VI]. A more thorough 
analysis of this line complex at slightly higher resolution does not find conclusive evidence for the 
[Fe II] line in this galaxy (Murphy et al. 2000a). 

Searching for the 1.9629 /im [Si VI] line was a primary motivation for the present spectroscopic 
survey. With a 167 eV excitation potential, this line indicates a high probability of the presence 
of an active nucleus. While it is possible to produce Si 5+ atoms via powerful shocks, Marconi 
et al. (1994) find that photoionization is the more likely excitation mechanism among a sample of 
nearby Seyfert galaxies. As stars are incapable of producing significant quantities of photons at this 
high energy, the natural photoionization source is an active nucleus with its power-law spectrum 
extending into the far-ultraviolet wavelength range. 

Only one galaxy in the sample of 33 has prominent [Si VI] emission. IRAS 08311—2459 has 
a very strong [Si VI] line, comparable in strength to the adjacent H~2 1-0 S(3) line, and appearing 
broader than the other lines in the spectrum. A detailed discussion of the [Si VI] line in this galaxy 
can be found in Murphy et al. (2000a) and Murphy et al. (1999). 

The only other galaxy in the sample with possible [Si VI] emission is IRAS 05246+0103. The 
spectrum of this galaxy shows emission at the expected position of [Si VI], with a signal-to-noise 
ratio of ~5. The adjacent H2 line appears broad on the short-wavelength side as well, calling into 
question the nature of the emission on the long-wavelength side of this line. No emission appears at 
this position in the secondary nucleus spectrum of Figure 2, ruling out a problem with atmospheric 
calibration or related issues. The H2 1-0 S(l) line appears unresolved in the primary nucleus, 
strongly suggesting that the [Si VI] emission is real. While IRAS 05246+0103 has no evidence of 
broad line emission, it does have a peculiar radio spectrum. Crawford et al. (1996) find this galaxy 
to be a gigahertz-peaked-spectrum (GPS) radio source, with a 6 cm flux density much higher than 
that measured at 20 cm. This condition is not found in many ULIRGs, with the Parks (PKS) galaxy 
IRAS 13451+1232 being the only other 2 Jy ULIRG identified as a GPS source. This galaxy may 
also have [Si VI] emission, as reported by Veilleux, Sanders, & Kim (1997), though the blend with 
the H2 line at their lower resolution makes it difficult to estimate the relative contribution of [Si VI] . 
The rest of the H2 lines in IRAS 13451+1232, as presented in Veilleux, Sanders, & Kim, indicate 
a H2 rotation temperature around 2500 K, in which case the [Si VI] line would account for ~30% 
of the total combined line flux. IRAS 13451+1232 does show impressive broadening of the Paa 
line, clearly standing out as an AGN-dominated source. While IRAS 05246+0103 has a similar 
radio spectrum to IRAS 13451+1232, it is almost three orders of magnitude less luminous at 20 
cm, according to Crawford et al. Perhaps there is a low-luminosity active nucleus component to 
IRAS 05246+0103, with the total infrared luminosity dominated by star formation. 
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4.4. Relative Nuclear Velocities 

The double nucleus ULIRGs often exhibit strong Paa emission from both nuclei, which can be 
used to study orbital dynamics of the galactic pair. Table 6 lists the relative radial velocities for 
the 12 ULIRGs with measured Paa velocities on both nuclei. The reported errors are considerably 
smaller than would be obtained by simply subtracting the velocities presented in Table 4 from 
those in Table 3, reflecting the fact that the relative velocities are not affected by systematic errors 
associated with wavelength calibration or slit illumination. 

Surprisingly, the velocity differences are very small, with the largest difference at 120 km s _1 . 
Random projection effects diminish the observed relative velocities statistically by a factor of two, 
as the most probable observing configuration involves pure angular motion on the plane of the 
sky. If all ULIRGs have the same relative nuclear velocity, random orientations would produce 
a distribution of observed velocities described by a cosine curve, with the peak at zero observed 
relative velocity. Binning the data from Table 6 into 50 km s _1 bins, the best fit to a cosine function 
has a maximum velocity of about 190 km s _1 , shown in Figure 8. Using 40 km s _1 bins results in 
a 165 km s _1 maximum relative velocity. 

The lack of high relative velocities among ULIRGs contrasts with the predictions obtained 
either from simple Newtonian calculations or from more sophisticated TV-body models of galaxy 
interactions which include gas dynamics and dynamical friction. Below we discuss these two ap- 
proaches, with attention to how these models might be modified to match the observations. 

First, we will treat the galaxies as point masses approaching each other from far away. These 
calculations assume conservative dynamics — i.e., that all loss of potential energy is converted into 
kinetic energy of the galaxies. In reality some fraction of the energy is lost to dynamical friction, 
which is why the galaxies do eventually merge. The foregoing discussion will not explicitly treat 
the frictional component, but will concentrate on the more basic elements of merger dynamics that 
offer much greater latitude in the ultimate relative velocities calculated. In the point mass scenario, 
the galaxies would reach a relative velocity Av = 2yf GM/r in excess of 1300 km s _1 at a separation 
distance of 10 kpc, assuming total system masses of 10 12 M Q per galaxy. Even ten times less mass 
per galaxy results in over 400 km s _1 of relative motion at 10 kpc. Changing the initial separation 
from infinity to 100 kpc or 50 kpc reduces these speeds by 5% and 10%, respectively. 

The high velocities computed above are obviously not consistent with the observations, but 
perhaps understandably given that galaxies are not well represented by point masses at close range. 
Galaxy masses are dominated by very large dark halos that act to soften the potential gradients — 
i.e., accelerations — experienced during a close encounter. 

Dark halos are thought to exist based on the anomalously flat rotation curves of galaxies to 
very large radii. Such phenomena suggest the enclosed spherical mass is proportional to radius, such 
that the circular rotation velocity, v c = \J GM enc i / r is constant. Rotation curves of large samples of 
galaxies via optical emission lines (e.g., Dale et al. 1999) indicate that this scaling persists at least 
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Fig. 7. — Rotation temperatures determined for the 33 primary ULIRG nuclei. The majority 
are consistent with thermal emission, for which the rotation temperature represents the ambient 
equilibrium temperature. A few are definitely mixtures of thermal and fluorescent processes, with 
vibration temperatures in the range of 3100-5000 K. 



Table 6. Relative Nuclear Radial Velocities 



Galaxy 


Relative Velocity 


Projected Separation 




(km s" 1 ) 


(kpc) 


IRAS 01521+5224 


-50 ± 25 a 


7.6 


IRAS 05246+0103 


-1 ± 12 


9.8 


IRAS 08572+3915 


-26 + 11 


5.8 


IRAS 09061-1248 


12 + 28 


7.2 


IRAS 09583+4714 


-53 ± 30 


37.5 


IRAS 10035+4852 


8 + 28 


11.2 


IRAS 10190+1322 


114 + 11 


5.3 


IRAS 12112+0305 


-98 + 13 


3.7 


IRAS 14348-1447 


-112 ± 10 


4.7 


IRAS 15245+1019 


-54 + 27 


3.4 


IRAS 17028+5817 


120 + 22 


23.0 


IRAS 22491-1808 


-3+10 


2.2 




a From Murphy et al. (2000b). The current dataset yields -22 + 
169. 
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out to 10-15 kpc. Observations in H I find a continuation of the fiat rotation curves at radii 2-3 
times the optical radii, often to 20-30 kpc (e.g., Casertano & van Albada 1990; Fich & Tremaine 
1991). Assuming a spherical halo distribution, the total halo mass following the M enc \ oc r scaling 
is given by 

M -- 10l '(208l^) 2 (^) Ms ' 

where v c is the circular velocity, and r cuto fr is the outer dimension of the halo distribution. Under 
this prescription, two galaxies falling from infinity are moving with a relative velocity of y/2v c by the 
time the halos touch, or when the galaxies are 2r cu t ff apart. The galaxies will clearly accelerate 
even more as the separation decreases to ~ 10 kpc, even though the acceleration is somewhat 
moderated by the effect of overlapping halo masses. This already presents a problem in that the 
observed relative velocities are less than the typical rotation velocities observed in ULIRGs. 

The simulations by Mihos & Hernquist (1996) provide a means of evaluating encounter ve- 
locities via models employing realistic halo potentials and dynamical friction. Output from these 
models is shown in Figure 9. The model defines a halo that yields a total galaxy mass of 4 x 10 11 
Mq when scaled to the Milky Way's rotation velocity. The halo has a relatively small physical 
extent, containing most of its mass within 25 kpc, and with an effective cutoff radius at 40 kpc. 
The simulations start with identical galaxies approaching from infinity (parabolic orbits), resulting 
in time-averaged velocities of ~ 350 km s _1 when the nuclei are between 4-25 kpc apart. This 
average sums over first, second, and third encounters. The peak velocities are in the neighborhood 
of 500 km s _1 , as seen in Figure 9. Starting the galaxies from rest at a 100 kpc separation rather 
than infinity reduces the peak velocities by only ~ 15%. 

Many ULIRGs are found at times well after the first major encounter, judging by the fact that 
their tidal features, formed at the epoch of the first significant encounter, are often vastly more 
extended than the nuclei are separated. However, some ULIRGs are believed to be experiencing 
their first major encounter (Murphy et al. 2000b). Looking at the first-encounter and late-encounter 
model predictions separately, the simulation discussed above predicts first-encounter velocities av- 
eraging 475 km s -1 , after which the galaxies separate by 50 kpc, followed by a final merger sequence 
with an average velocity around 275 km s _1 , ranging from 0-500 km s" 1 during this time. 

Models with larger-scale halos maintaining the same internal rotation structure demand higher 
total masses, resulting in increased net velocity differences between the merging galaxies. On the 
other hand, if ULIRGs lacked any halo material past 10 kpc, corresponding to total masses around 
10 11 M Q , the galaxies would behave much more like point masses, and dynamical friction would be 
less important. In this scenario, the 10 kpc velocities would be approximately 400 km s _1 . 

It is very difficult to reproduce the small observed relative velocities among ULIRGs with simple 
physical models — with or without large halos. ULIRGs are massive, gas rich galaxies with rotation 
velocities and near-infrared absolute magnitudes characteristic of galaxies at least as massive as 
the Milky Way. The highest observed relative velocity in this sample is a mere 120 km s _1 , 
though the observed distribution is consistent with a typical deprojected velocity difference of 200 
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Fig. 8. — Histogram of the absolute values of projected relative nuclear velocities. If all ULIRGs 
had the same relative velocity, projection effects would result in our observing a cosine distribution. 
The best fit cosine to the data has a 190 km s~ x velocity amplitude. 
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Fig. 9. — Simulations of galaxy interactions, courtesy of J. C. Mihos. At left is the orbital track, at 
center is nuclear separation as a function of time, and at right is the relative velocity as a function 
of time. The model galaxies each have a total mass of 4 x 10 11 M©, and possess Milky- Way-like 
disk rotation velocities. 
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km s 1 . Some ULIRGs have been observed with velocity differences very near this value (e.g., 
IRAS 20046-0623; IRAS 19252-7245: Murphy et al. 2000b; Mihos & Bothun 1998), but none are 
seen near the 400 km s _1 value as expected from a variety of physical predictions. 

One possibility deserves mention, which is that during the late-merger state, eccentric orbits 
may conspire with the delay in coordinating the accumulation of fuel for a nuclear starburst in such 
a way as to select ULIRGs at the apocenters of their orbits, when velocities are low. Referring to 
Figure 9, during the final merger process (well after the initial encounter), the nuclear orbits are 
very eccentric, owing to the removal of angular momentum, largely by tidal debris. The very high 
velocity peaks near 500 km s _1 are associated with the close passages. If the ultraluminous activity 
is episodic in nature, with relatively short duration bursts occurring in reaction to transitory tidal 
disruptions, then the ULIRGs we see may represent stages shortly after these close encounters. 
A finite time is required to organize gas motions in response to the close encounter, with the 
ultraluminous activity turning on only after a significant gas concentration has been established. 
It is not difficult to imagine in this scenario that the gas — moving more slowly than the galaxy, 
but with less distance to travel — will reach the nucleus at approximately the same time that the 
galaxies reach apocenter. Keeping in mind that we are discussing the events occurring after the 
second close encounter, the orbital timescale is ~ 10 8 yr, such that apocenter is reached a fewxlO 7 
yr after close approach. This timescale is in rough agreement with expectations of the dynamical 
timescale within a few kpc of the nucleus — the region from which the fuel for the nuclear event is 
gathered. During the final merger process, it is possible that several ultraluminous bursts occur 
following close interactions of the nuclei. The idea of repetitive starbursts in the final stages of 
merging was first explored by Noguchi (1991), though in the context of cloud-cloud collisions in this 
case. The scheme proposed here of delayed — and potentially episodic — bursting would bias late- 
merger ULIRGs to larger separations and therefore lower velocities, quite possibly consistent with 
the observed distribution. The problem with first encounter ULIRGs remains, since the average 
first encounter velocity is around 475 km s _1 . Unfortunately the current sample is not large enough 
to probe possible differences in the velocity distributions of early and late ULIRG encounters. 



5. Summary 

Near-infrared spectra in the 2 ^m window have been obtained for a nearly complete, volume- 
limited sample of 33 ULIRGs. Lines of atomic hydrogen and helium recombination, molecular 
hydrogen vibration-rotation, [Fe II], and [Si VI] are seen in the spectra. The Pact is predominantly 
seen as a nuclear event, with the majority of the emission arising in the central kpc of most ULIRGs. 
Initiated as a search for optically hidden active galactic nuclei, this survey reveals only two such 
galaxies, both of which are characterized as Seyfert galaxies based on visible-light spectra. It is 
found via comparison of the Paa and Br7 lines that significant extinction of several magnitudes 
persists into the near-infrared bands, such that buried active nuclei can not be ruled out in many 
objects. It is estimated that with the sensitivity of the current survey, at least 25 mag of visual 
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extinction would be needed to hide a QSO with half the bolometric luminosity of the ULIRG. 

Strong molecular hydrogen emission is seen in all of the primary nuclear spectra. The molecular 
hydrogen spectrum indicates that most H 2 emission in ULIRGs stems from thermal excitation 
rather than fluorescent processes. Some ULIRGs appear to have a small fraction of their molecular 
hydrogen emission owing to fluorescence, though none approach pure fluorescence. 

The relative radial velocities found in this sample of ULIRGs does not exceed ~ 120 km s _1 , 
with an inferred maximum deprojected velocity differential of around 200 km s _1 . This relatively 
low velocity differential is difficult to produce from interactions between galaxies with masses com- 
parable to that of the Milky Way using simple physical assumptions. It is suggested that ULIRGs 
in the late stages of merging may be selectively found near the time of apocenter, leading to a lower 
average relative velocity between nuclei. 
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A. Procedures for Measuring Line Properties 

The line properties reported in Tables 3 and 4 represent measured values from the spectra, 
rather than Gaussian fit parameters. Especially for the weaker lines, the quality of the continuum 
subtraction is vital to proper line characterization. This is done with a simultaneous multiple 
Gaussian line fit and quadratic continuum fit. In all but a few cases, this results in a very close fit 
to the continuum. For those few cases where the continuum fit is questionable, some of the weaker 
lines are excluded from Tables 3 and 4. All of the fitting procedures make use of the noise estimate 
computed on a pixel-by-pixel basis, as described earlier. In the end, the noise is scaled slightly 
(typically by factors of less than 30%) in order to make the final reduced chi-squared value equal 
to one. The scaled noise is used in estimating errors on the values in Tables 3 and 4. Each of the 
line parameters is computed based on the following line definition procedure. Each potential line is 
searched for by summing two resolution elements (7 pixels) of information around the line center, 
and determining if there is emission above the continuum exceeding 2.5a. For lines satisfying this 
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condition, cutoffs are established on the red and blue sides of the line center by the following rules: 

1. If the spectrum extends below the continuum level, establish a cutoff, including only positive 
valued points. 

2. If there exists a local minimum, followed by at least three higher points, assume that the 
following points belong to another line and exclude them, keeping only the minimum point. 

3. If the next two pixels are between 0-la above the continuum, include them and quit. 

These procedures obviously can not practice the same discrimination that a diligent scientist might, 
but we find that a reasonable job is performed — especially on lines whose boundaries are obvious 
to the eye. Nonetheless, lines with a final computed equivalent width measuring less than 2.5a 
are typically excluded, though those that appear to be real in the one-dimensional spectra are 
preserved. Conversely, a few lines with reported detections greater than 2.5a are excluded on the 
basis of either poor local continuum fit or simply less-than-believable results. 

Because the weaker lines do not have a well defined full-width at half-maximum, the FWHM 
values in Tables 3 and 4 are based on a combination of line flux and peak amplitude. If a Gaussian 
shape is assumed, the FWHM of the line would be equal to the line flux divided by the peak 
amplitude, times the proportionality constant of ^/ln(256) /2ir = 0.94. Performed in this manner, 
it is relatively simple to compute the error of such an estimate. A strict measure of the FWHM is 
provided in Tables 3 and 4 for the Paa line in each galaxy for comparison to the value estimated 
by this technique. All of the reported FWHM values have been reduced to intrinsic line widths 
by subtracting the instrumental resolution in quadrature. Each spectrum uses the measured line 
widths in the wavelength calibration spectrum as a basis for the decomposition. Values in Tables 3 
and 4 for which the measured line width is exceeded by the instrumental width are listed as 
unresolved lines of zero width. Almost all of the lines in this survey are resolved, typically around 
200-300 km s" 1 in width. 

The velocity measurements in Tables 3 and 4 are indicated relative to Paa, though the Paa 
velocity itself is expressed as a direct quantity in km s _1 . These values are simply computed as 
first-moment measures in the spectral lines, with error estimates computed from the noise estimates 
in a straightforward way. Added to this error estimate is a wavelength calibration error, assumed 
to be 10-15 km s _1 , and a slit placement error. The latter uncertainty attempts to estimate 
the wavelength offset induced by a non-uniform slit illumination. When the line source is spatially 
compact enough to produce a discernible peak somewhere in the slit, the location of this peak within 
the slit influences the ultimate wavelength associated with this emission. For these purposes, it is 
assumed that the line source has a real physical size roughly half that of the continuum, which is 
almost always spatially resolved. Further, it is assumed that the source is centered in the slit with a 
Gaussian la placement error of one-quarter of a slit width. Only cases of exceptional seeing result 
in an appreciable error from this source, with a median error estimate of 15 km s _1 and a maximum 



-47- 



of 38 km s . Despite the accounting of uncertainties, the velocities indicate a dispersion outside 
of the error budgets that is not correlated with line species or grating setting. It may well be that 
the variable spectral range over which the line is summed is largely responsible for determining the 
offset, as the range is biased to contain only positive points, and will extend as far as possible from 
the line center to include such points. With the large moments generated from distant points, the 
estimated velocity may wander far from its true center. Thus the values of velocity in Tables 3 and 
4, especially for weaker lines, are perhaps of little use. Gaussian fit centers may do a better job of 
appropriately representing the line velocities. 



B. Estimating Extinction to Hidden AGN 

The median QSO spectrum is constructed from nine optically identified QSOs in the redshift 
range 0.089 < z < 0.182, seven of which have available K band photometry. Comparing the 
absolute K band magnitudes to bolometric luminosities of the quasars, one finds the approximate 
relation Mr = 4.0 ± 0.5 — 2.5 log (L^/Lq). The median infrared luminosity of the 31 non-AGN 
ULIRGs in the present sample is 1.1 x 10 12 L Q . A QSO within a ULIRG that contributes half 
of the bolometric luminosity — most of which is reprocessed into the infrared — would then have 
Lbol = 5.5 x 10 11 Lq, or Mk = —25.4 ± 0.5 mag. This establishes the K band flux density of the 
QSO continuum. At a rest wavelength of 1.881 ^m, corresponding to +1000 km s _1 in velocity, it 
is found that the median QSO has an excess flux density due to Paa measuring 22% of the baseline 
continuum level (Murphy et al. 1999). Thus the absolute magnitude of the line emission at this 
velocity is M K (Paai o ) = -23.74 ± 0.5 mag. 

Typical ULIRG continuum levels can be approximated in a similar way as that above, using 5" 
aperture photometry of the sample galaxies where available (Kim 1995; Murphy et al. 2000b). Eight 
of the 31 non-AGN ULIRGs have K s measurements, and comparing to the infrared luminosity (8— 
1000 ^m) leads to the relation Mk s = 6.3 ±0.6 — 2.5 \og{Lm/ Lq). Therefore the typical non-AGN 
ULIRG in this sample, with a median luminosity log(L/ij/L ) = 12.05, has Mk s = —23.83 ± 0.6 
mag. 

The expected ratio of broad Paa to the ULIRG continuum at 1.881 ^m is then 0.92, with the 
la range extending from 0.44 to 1.91. If a given ULIRG has a spectrum with a signal-to-noise ratio 
G, then diminishing the QSO Paa to the la level requires 

Qpg a e~ T _ J_ 

Qcont^ T ~\~ U con t G 

where Q represents QSO contributions to continuum and line flux, U is the ULIRG continuum 
contribution, and e~ T is the attenuating factor. With the additional relations that Qpaa/Qcont = 
0.22 and Qp a a/U con t = 0.92, we find that the required attenuation is e T = 0.92(G — 1/0.22), or 
2.51og(0.92G — 4.2) magnitudes at Paa. The visual extinction is roughly 7 times this value. 
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C. Molecular Hydrogen Diagnostics 

C.l. Population Diagrams 

The specific mechanism responsible for the H2 emission can be investigated via population 
diagrams. By accounting for various degeneracies, transition probabilities, and photon energy, one 
may compute the relative populations of the upper states of each transition based on the observed 
line flux. In pure thermal equilibrium, the population levels will follow an exponential profile with 
respect to upper state energy. More formally, the intensity of a given transition is described by 

-Ex 

I(v 1 ,v 2 ,Ji,J2) cx g s {Ji)A(v 1 ,v 2 , Ji, J 2 ){2Ji + l)hu 12 e w, 

where the upper and lower states are denoted by the 1 and 2 subscripts, v represents vibration state, 
J the rotation state, g s {J\) the ortho/para degeneracy, A(v±, v 2 , Ji, J 2 ) the transition probability, 
v\ 2 the emitted photon frequency, E\ the upper state energy, and T the equilibrium temperature. 
Therefore a logarithmic plot of the line intensities divided by the terms multiplying the exponential 
above produces a linear relation against the upper energy scale with a slope proportional to the 
inverse of the temperature. The physical parameters for selected K band transitions are given in 
Table 7. 

When the molecular hydrogen is not in thermal equilibrium — such as may be the case in 
fluorescent excitation by ultraviolet photons — a diagram such as that described above will not 
show all points lying along a common line, but rather disconnected linear segments, with common 
vibration states each defining a line connecting the different rotation states. The characteristic 
temperature of this line is referred to as the rotation temperature, T ro t. Ideally, each vibration state 
has the same characteristic rotation temperature, and these segments can be connected with another 
line characterized by a vibration temperature, T v n,. The thermal case simply has T ro t = T v n, = T. 
Typically, pure fluorescent conditions are characterized by T rot = 800 — 1400 K, and T vi b = 5000- 
6000 K (Black k, van Dishoeck 1987; Takayanagi, Sakimoto, & Onda 1987). Excitation conditions 
associated with thermal processes are typically found to have T rot = 1600-2400 K, with a slight 
difference between rotation and vibration temperatures characterized by T ro t ~ 0.81^, (Tanaka et 
al. 1989; Hasegawa et al. 1987; Beckwith et al. 1983). 

C.2. Typical K Band H 2 Line Ratios 

Because ULIRGs tend to show relatively uniform H2 emission properties — namely those cor- 
responding to thermal emission at 1500-2500 K — it may be useful for other researchers to have a 
table of expected H2 line strengths in the K band for a range of typical temperatures. Table 8 
provides such a reference, with line intensities normalized to the H 2 1-0 S(l) line. These estimates 
assume an ortho-to-para ratio of 3:1, and T ro t = 0.81^,, as is found empirically (Tanaka et al. 1989; 
Hasegawa et al. 1987). In order to regain the pure thermal regime characterized by a single tern- 
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Table 7. Physical Parameters of H2 Transitions 



Transition 




Ji 


A(vi,V2, Ji, J 2 ) a 


£i/fc b 




w 




(10- 7 s- 1 ) 


(K) 


1-0 S(0) 


2.22268 


2 


2.53 


6471 


1-0 S(l) 


2.12125 


3 


3.47 


6951 


1-0 S(2) 


2.03320 


4 


3.98 


7584 


1-0 S(3) 


1.95702 


5 


4.21 


8365 


1-0 S(4) 


1.8914 


6 


4.19 


9286 


1-0 S(5) 


1.8353 


7 


3.96 


10341 


2-1 S(0) 


2.3550 


2 


3.68 


12095 


2-1 S(l) 


2.2471 


3 


4.98 


12550 


2-1 S(2) 


2.1536 


4 


5.60 


13150 


2-1 S(3) 


2.0729 


5 


5.77 


13890 


2-1 S(4) 


2.0035 


6 


5.57 


14764 


2-1 S(5) 


1.9443 


7 


5.05 


15763 


2-1 S(6) 


1.8942 


8 


4.30 


16880 


2-1 S(7) 


1.8523 


9 


3.38 


18107 


3-2 S(3) 


2.2008 


5 


5.63 


19085 



a Taken from Turner, Kirby-Docken, & Dalgarno. (1977). 
b Taken from Dabrowski (1984). 
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perature, the 2-1 entries in Table 8 should be multiplied by exp(— 1107/T), and the 3-1 transition 
multiplied by exp(— 2144/T), lowering these values considerably 
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Table 8. Expected H 2 Line Strengths Relative to H 2 1-0 S(l) 



Transition \ air (/an) T rot = 1500 K T rot = 2000 K T rot = 2500 K T rot = 4000 K 



1-0 S(0) 


2.22268 


1-0 S(l) 


2.12125 


1-0 S(2) 


2.03320 


1-0 S(3) 


1.95702 


1-0 S(4) 


1.8914 


1-0 S(5) 


1.8353 


2-1 S(0) 


2.3550 


2-1 S(l) 


2.2471 


2-1 S(2) 


2.1536 


2-1 S(3) 


2.0729 


2-1 S(4) 


2.0035 


2-1 S(5) 


1.9443 


2-1 S(6) 


1.8942 


2-1 S(7) 


1.8523 


3-2 S(3) 


2.2008 



0.228 


0.211 


1.000 


1.000 


0.336 


0.374 


0.805 


1.019 


0.177 


0.261 


0.295 


0.519 


0.015 


0.030 


0.068 


0.143 


0.023 


0.053 


0.055 


0.145 


0.012 


0.037 


0.020 


0.072 


0.003 


0.014 


0.004 


0.020 


0.003 


0.016 



0.201 


0.187 


1.000 


1.000 


0.398 


0.438 


1.174 


1.450 


0.329 


0.468 


0.728 


1.211 


0.045 


0.083 


0.225 


0.441 


0.089 


0.191 


0.259 


0.622 


0.072 


0.197 


0.156 


0.496 


0.033 


0.124 


0.054 


0.246 


0.045 


0.202 



Note. — This table assumes that T ro t = 0.8T v n,. To recover a strictly thermal regime (T ro t = 
T vib = T), multiply the 2-1 transistions by exp(-1107/T) and the 3-2 transition by exp(-2144/T). 
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Table 3. Measured Line Properties for Primary Nuclei 



Galaxy 


Continuum 


Pan FWHM b 


Rows c 


He I 


Pact 


H 2 S(3) 


[Si VI] 


[Fe II] 


H 2 S(l) 


Br7 




Slope a 


(km s" 1 ) 




1.8689 


1.8751 


1.9570 


1.9629 


1.9670 


2.1213 


2.1655 



IRAS 00262+4251 



IRAS 01521+5224S -1.94 



IRAS 04232+1436 



IRAS 05246+0103B 1.18 



IRAS 08030+5243 



IRAS 08311-2459 



IRAS 08344+5105 



IRAS 08572+3915N 3.26 



IRAS 09061-1248S -3.28 



IRAS 09111-1007 



IRAS 09583+4714NB -3.16 



IRAS 10035+4852NB -4.78 



BW 
FWHM 



BW 
FWHM 



BW 
FWHM 



BW 
FWHM 



BW 
FWHM 



BW 
FWHM 



BW 
FWHM 



BW 
FWHM 



BW 
FWHM 



BW 
FWHM 



BW 
FWHM 



BW 
FWHM 



0.42+0.1 

420+240 
-67+78 



1.09+0.13 

497+78 
-83+51 

0.53+0.13 
340+160 
97+71 



0.67+0.13 
60+110 
-45+39 

0.84+0.13 
440+110 
-51 + 53 

1.16+0.5 

240+280 
-260+120 

0.42+0.11 
290+210 
67+99 

0.74+0.11 

0+63 
-78+24 

0.454+0.049 
0+61 

-43+29 



3.52+0.15 

366+27 
29204+18 

12.02+0.44 

356+23 
24011 + 16 

8.74+0.25 

218+18 
23666+15 

8.89+0.19 
190.1+9.6 

29105+22 

14.46+0.16 
326.2+6.8 
24959+12 

22.69+0.29 

477+10 
30114+16 

20.08+0.26 

300+13 
29002+12 

11.192+0.1 
334.8+8 
17493+19 

10.33+0.25 

320+24 
22017+13 

5.53+0.1 
246+13 
16231 + 12 

16.75+0.2 
263.6+6.1 
25754+12 

8.71+0.19 

254+12 
19435+17 



2.95+0.14 
480+42 
-18+29 

0.79+0.14 
170+96 
-58+46 

1.34+0.15 
195+70 
25+33 

3.15+0.26 

575+77 
255+47 

1.3+0.11 

331 + 77 
-29+34 

3.68+0.15 
501 + 41 
129+26 

1.64+0.13 
196+40 
60+28 

0.695+0.062 
139+45 
-69+37 

1.215+0.092 
421 + 99 
-81 + 26 

0.736+0.033 
189+27 
-12+17 

3.58+0.22 
411 + 74 

25+25 

2.46+0.13 
315+43 
38+35 



0.404+0.08! 
260+160 
187+63 



0.86+0.2 
310+260 
-180+55 



3.005+0.088 
708+52 
78+20 



0.39+0.1 

0+91 
196+87 

0.488+0.078 
360+130 
337+45 



0.56+0.19 
100+240 
-72+96 

0.295+0.052 
310+91 
244+51 



0.21 + 0.11 

0+240 
180+170 



3.88+0.24 
380+56 
-44+27 

2.44+0.36 
290+120 
-101+34 

1.76+0.16 
203+51 
-11 + 26 

1.53+0.14 
161+40 
27+38 

1.36+0.22 
270+130 
-47+31 

2.734+0.092 
351 + 27 
104+21 

2.16+0.27 
310+100 
0+30 

0.448+0.027 
135+30 
-78+31 

1.53+0.13 

266+58 
34+23 

0.682+0.035 
164+25 
78+20 

4.06+0.25 
384+59 
13+28 

2.915+0.089 
256+18 
-2+28 



1.38+0.49 
0+110 
-95+76 

1.4+0.26 
220+100 
63+46 

0.61 + 0.24 
210+210 
70+110 

2.32+0.28 
380+110 
-173+44 

2.18+0.25 
470+170 
94+34 

2.32+0.36 
310+140 

52+35 

0.762+0.072 
271+62 
-56+37 

1.47+0.23 
260+110 
7+45 

0.59+0.062 
168+66 
54+26 

1.85+0.25 

256+98 
-134+38 

1.08+0.15 
201 + 71 
46+44 



Table 3 — Continued 



Galaxy 


Continuum 


Pact FWHM b 


Rows c 


He I 


Paa 


H 2 S(3) 


[Si VI] 


[Fe II] 


H 2 S(l) 


Br7 




Slope a 


(km s" 1 ) 




1.8689 


1.8751 


1.9570 


1.9629 


1.9670 


2.1213 


2.1655 



IRAS 10190+1322B -1.83 



IRAS 10494+4424 -3.98 



IRAS 11095-0238 



IRAS 12112+0305SW 1.17 



IRAS 14348-1447SW -0.98 



IRAS 14352-1954 



IRAS 14394+5332 



IRAS 15245+1019B -2.34 



IRAS 15250+3609 



IRAS 15462-0450 



IRAS 16487+5447W 1.47 



IRAS 17028+5817W -0.61 



IRAS 18470+3233 



BW 
FWHM 



BW 
FWHM 



BW 
FWHM 



BW 
FWHM 



BW 
FWHM 



BW 
FWHM 



BW 
FWHM 



BW 
FWHM 



BW 
FWHM 



BW 
FWHM 



BW 
FWHM 



BW 
FWHM 



BW 
FWHM 



1.04+0.11 
510+160 
-51+34 

1.62+0.3 
360+180 
-38+64 

1.68+0.61 
380+340 
80+140 

1.47+0.15 

550+120 
150+43 

1.102+0.089 

258+55 
19+33 

1.602+0.064 
558+77 
-137+33 

0.277+0.069 
0+75 
25+39 

0.43+0.23 

0+180 
-70+200 



1.96+0.18 
500+140 
29+40 

0.631+0.097 
320+130 
46+46 

2.21+0.4 
200+120 
-9+51 



9.6+0.13 
437+16 
22889+19 

14.96+0.12 
404.2+7.1 

27675+16 

20.38+0.39 

286+12 
31968+18 

13.97+0.39 

303+21 
21980+21 

12.78+0.19 
236.8+8.4 

24802+22 

23.54+0.23 
205.3+4.9 
26942+17 

15.75+0.13 
381.8+8.7 
31373+21 

15.32+0.14 
294.6+7.3 
22634+16 

9.49+0.16 
147.8+9.9 
16535+17 

5.232+0.085 
535+15 
29917+24 

29.02+0.34 
332.3+9.1 
31106+23 

11.89+0.13 
315.6+8.7 
31805+24 

41.07+0.44 
258.2+6.9 
23517+21 



1.164+0.081 
412+72 
12+37 

1.612+0.072 
377+38 
-2+29 

5.65+0.25 
350+49 
-59+26 

1.095+0.086 

295+63 
-65+42 

3.15+0.12 

450+35 
-65+35 

1.19+0.11 

193+58 
9+37 

5.03+0.097 

490+22 
-5+31 

1.324+0.067 
307+46 
-19+30 

1.45+0.14 

323+88 
6+35 

0.48+0.029 

297+53 
229+40 

7.13+0.29 
374+36 
-52+33 

1.26+0.078 

250+48 
3+33 

2.04+0.16 

305+68 
11+44 



0.154+0.072 
0+130 
-50+140 



0.34+0.14 

220+290 
-100+190 

0.244+0.083 
120+210 
170+110 



0.143+0.08 

160+280 
210+170 

0.341 + 0.078 
300+210 
137+67 

0.55+0.1 

400+180 
-15+60 



0.35+0.13 
140+170 
120+130 

0.41+0.16 

200+300 
140+120 



1.248+0.096 

336+65 
7+32 

1.89+0.11 

427+65 
-3+29 

8.56+0.61 
241+43 
-73+35 

0.957+0.071 
0+31 
-19+30 

3.11 + 0.14 

334+36 
-46+30 

1.88+0.15 

243+48 
35+34 

4.9+0.16 

426+33 
61 + 33 

1.242+0.088 
247+50 
25+27 

1.928+0.059 
223+19 
58+32 

0.268+0.024 
102+43 
158+32 

5.95+0.39 
265+54 
-134+30 

1.53+0.14 
301 + 74 
6+39 

2.46+0.17 
319+49 
11 + 34 



1.12+0.11 

320+91 
22+34 

1.73+0.16 

364+75 
-51+31 

2.5+1 
0+100 
80+130 

1.57+0.15 

185+53 
-51+33 

1.35+0.2 

127+80 
35+40 

2.97+0.26 
166+46 
6+32 

2.14+0.27 
370+110 
127+47 

2.08+0.13 
304+51 
-15+25 

1.429+0.077 
204+34 
100+33 

0.516+0.06 
360+110 
237+44 

3.46+0.89 
230+140 
-195+78 

2.13+0.25 

336+94 
22+45 

3.52+0.25 
249+50 
-101 + 30 



Table 3 — Continued 



Galaxy 


Continuum 


Pact 


FWHM b 


Rows c 


He I 


Paa 


H 2 S(3) 


[Si VI] 


[Fe II] H 2 S(l) 


Br7 




Slope a 


(km s" 1 ) 




1.8689 


1.8751 


1.9570 


1.9629 


1.9670 2.1213 


2.1655 


IRAS 19458+0944 


-0.25 




269 
























BW 


0.668+0.086 


11.093+0.075 


1.46+0.13 




1.236+0.089 


1.74+0.29 










FWHM 


390+180 


272.9+5.1 


324+84 




237+50 


150+110 










cz 


-27+58 


29982+35 


-83+59 




3+47 


73+53 


IRAS 20046-0623 


-3.99 




244 


BW 
FWHM 

cz 




6.3+0.21 

247+27 
25293+14 


1.16+0.2 
190+100 
-16+43 




0.45+0.2 1.63+0.44 
0+140 0+88 
50+140 6+60 


3.67+1 
310+170 
100+80 


IRAS 20087-0308 


-0.3 




595 
























BW 


0.465+0.093 


10.22+0.14 


1.611 + 0.067 




2.11 + 0.15 


0.64+0.24 










FWHM 


270+150 


574+23 


475+49 




570+110 


280+280 










cz 


-92+48 


31600+20 


-41 + 33 




64+30 


36+97 


IRAS 20414-1651 


-1.93 




515 
























BW 


0.82+0.16 


10.47+0.21 


1.64+0.16 




1.71 + 0.17 


2.23+0.63 










FWHM 


400+190 


488+20 


287+82 




460+120 


170+180 










cz 


53+71 


26043+25 


16+49 




35+46 


-137+92 


IRAS 21504-0628 


-1.29 




256 
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FWHM 


230+100 


260.8+5.6 


250+130 




310+43 


144+40 










cz 


-17+52 


23352+26 


46+44 




43+43 


-4+41 


IRAS 22491-1808B 


-4.12 




128 


BW 
FWHM 

cz 




13.21 + 0.38 
131 + 18 

23277+24 


3.06+0.23 
0+52 
-35+41 




5.26+0.37 
302+60 
-2+33 


3.04+0.83 
430+350 
10+77 


IRAS 23327+2913S 


0.08 




151 


BW 
FWHM 

cz 


0.43+0.11 
170+130 
68+76 


6.85+0.12 

193+9 
32235+18 


0.87+0.26 
110+120 
-1 + 77 




1.92+0.36 
210+120 
219+51 




IRAS 23365+3604 


-1.46 




186 


BW 
FWHM 

cz 




9.284+0.083 
191.1 + 5 
19309+17 


1.291 + 0.079 
245+45 
60+24 




1.493+0.057 
219+22 
-36+25 


1.345+0.094 

216+39 
-9+28 



bo 
-4 



a Slope expressed as power-law index, a, where fx — C\ a . 

b Strictly formal measure of the full-width at half-maximum of the Paa line, deconvolved by the instrumental resolution. 

c For each galaxy, the following measures are provided for selected lines: equivalent width (EW) is given in nanometers; FWHM as estimated via line flux and peak 
amplitude (see text), deconvolved by the instrumental resolution; first moment centroid of the line, expressed as recessional velocity (cz) for Paa, and velocity offsets from 
Paa for all other lines. The first moment measure looses accuracy for weak lines depending on the spectral range summed, and therefore is of limited use here. 



